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Abstract 

Extending the statistical approach proposed in a parallel paper l|], 
propose a stochastic inverse kinetic theory for small-scale hydrodynamic turbulence based on the 
introduction of a suitable local phase-space probability density function (pdf ). In particular, we pose 
the problem of the construction of Fokker-Planck kinetic models of hydrodynamic turbulence. The 
approach here adopted is based on the so-called IKT approach (inverse kinetic theory), developed 
by Ellero et al. (2004-2008) which permits an exact phase-space description of incompressible fluids 
based on the adoption of a local pdf. We intend to show that for prescribed models of stochasticity 
the present approach permits to determine uniquely the time evolution of the stochastic fluid fields. 
The stochastic-averaged local pdf is shown to obey a kinetic equation which, although generally 
non-Markovian, locally in velocity-space can be approximated by means of a suitable Fokker-planck 
kinetic equation. As a side result, the same pdf is proven to have generally a non-Gaussian behavior. 

PACS numbers: 47.10.ad,47.27.-i,05.20.Dd 
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I. INTRODUCTION 

In this paper - extending the results of Ref.jlJ - we intend to formulate an IKT (in- 
verse kinetic theory) approach for the full set of fluid equations describing the phenomenon 
of turbulence in an incompressible fluid, here identified with the stochastic incompressible 
Navier-Stokes equations. Unlike Ref. where the set of stochastic-averaged fluid equa- 
tions were considered, we intend to show that the kinetic description, for the stochastic local 
probability distribution function (stochastic local pdf), is generally non-Markovian when the 
complete set of fluid equations is considered. The theory is shown to satisfy an H-theorem 
assuring the strict positivity of the local pdf as well of its stochastic average. In particular, 
the stochastic-averaged pdf is shown to satisfy an approximate Fokker-Planck-type kinetic 
equation. The explicit representation of the leading Fokker-Planck coefficients is provided. 
The result, which holds in principle for arbitrary prescribed stochasticity of the fluid fields, 
is achieved by means of an IKT which delivers the complete set of fluid fields, all expressed 
in terms of appropriate moments of the stochastic local pdf. 



A. Background and open problems 

The investigation hydrodynamic turbulence in incompressible fluids is nowadays playing 
a major role in fluid dynamics research. In fact the phenomenon of turbulence is essen- 
tially ubiquitous, being related to its statistical character. Indeed in many cases the fluid 
fields which define an incompressible isothermal fluid, i.e., the fluid velocity and pressure 
Z(r,t) ={\,p} are actually not known deterministically but only in a statistical sense. 
This means that the fluid fields must contain some kind of parameter-dependence Z(r, t, a), 
where a = (ai, a n ) G V a C M. n denotes a suitable stochastic real vector independent of 

(r,i) to which a stochastic probability density g(a) can be attached, so that J d n ag(a) = 1 

v a 

and furthermore the stochastic average (Z(r,t,a)) = J d n ag(a)Z(r,t,a) exists. As a 

v a 

consequence, the fluid fields can be represented in terms of the stochastic decomposition 
Z(r,t, a) = (Z(r,t)) + 5Z(r,t, a), where 5Z(r,t, a) = {<5V, 5p} are suitable stochastic fluctua- 
tions of the fluid fields. The precise form which these fluctuations may take defines what is 
usually denoted as the stochastic behavior (or stochasticity) of the fluid. The vector a, which 
spans a suitable subset V a of M™, can in principle be assumed either continuous or discrete. 
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Its possible definition, as well the identification of the related probability density appearing 
in the stochastic-averaging operator (), is however manifestly non-unique. In fact, these 
definitions are closely related to the types of stochastic behavior which may appear in the 
fluid. In the mathematical theory of turbulence one can distinguish in principle two possible 
types of stochasticity: either intrinsic or numerical. The first type (intrinsic stochasticity) 
arises when, leaving unchanged the functional form of the fluid equations, the fluid equations 
are intended as stochastic pde's. This happens if at least one of the following sources of 
stochasticity is introduced: 1) Stochastic initial conditions: in this case the initial fluid fields 
Z(r,t ) = Z Q (r) are assumed stochastic, i.e. of the form, Z Q (r,a) = (Z Q (r,a)) + 5Z Q (r,a), 
being (Z (r,a)) and 5Z (r,a) suitable vector fields. 2) Stochastic boundary conditions: this 
happens if the boundary fluid fields Z w (r,t)\ sn are prescribed in terms of a suitable stochas- 
tic vector field of the form Z w (r,t,a)\ sn = (Z w (r,t,a))\ sn + 5Z w (r,t,a)\ sa . Here, Dirichlet 
(no-slip) boundary conditions have been imposed (for the fluid fields) on the boundary 
set 5il of the fluid domain f2 C IR 3 by letting Z(r,t)\ sn = Z w (r,t)\ m . 3) Stochastic forc- 
ing. In this case the volume force density acting on the fluid is assumed stochastic, i.e., 
of the form f(r,t,a) = (f(r,t, a)) + 5f (r,t, a), being (f(r,t,a)) and 5f(r,t,a) suitable vector 
fields. It is obvious that the parameters a = (ai, a n ) and the related probability density 
g(a) can be set, in principle, arbitrarily. In fact, no information on them can be gath- 
ered from the deterministic fluid equations. For the same reason, also the definition of the 
stochastic fluctuations appearing in the previous equations, namely SZ Q (r, a), 5Z w (r,t, a)\ sn 
and Sf(r,t, a), remains unspecified. As a consequence, each of the stochastic vector fields 
,oz), Z w (r,t, ct) and f (r,t, o;) may in principle be characterized by different stochas- 
tic parameters a and probability densities g(a). Regarding, in particular, the definition 
of the boundary conditions, we remark that the stochasticity of the boundary fluid fields 
Z w (r,t,a)\ sn may be simply a result of the choice of the prescribed boundary 5Q. This 
happens if SQ is identified, for example, with a moving surface in which each point (of the 
surface) moves with random motion. The second type (numerical stochasticity) arises - in- 
stead - as a result of the approximate numerical solution methods adopted. All numerical 
methods, in fact, involve in some sense the introduction of appropriate approximations for 
the relevant differential operators, based on suitable time and space discretizations. As a 
consequence, the numerical solutions obtained for the fluid fields become inaccurate on scale 
lengths comparable or smaller than the spatial discretization (or grid) scales, thus producing 
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stochastic error fields. As a result, even if the fluid equations are treated as deterministic 
(which means ignoring all possible sources of intrinsic turbulence indicated above), nu- 
merical errors produce fluid fields which behave effectively as stochastic, giving rise to the 
phenomenon of numerical turbulence, also known as small-scale or sub-grid turbulence. This 
means that in principle it is possible to treat the two problems in a similar way. However, 
numerical stochasticity, as opposite to intrinsic stochasticity, is expected - in principle - to 
allow a well-defined statistical description. Nevertheless, a consistent theoretical formulation 
of small-scale turbulence of general validity seems still far away. In particular, still miss- 
ing is a consistent statistical description, based on the proper definition of the stochastic 
probability density g(a), able to recover from first principles the correct form of the pdf 
appropriate for arbitrary turbulence regimes. A widespread picture (of turbulence) consists 
both of an ensemble of finite-amplitude waves with random phase. However, there is an 
increasing evidence that this picture is an oversimplification. In fact, it is well known that 
turbulence may include fluctuations whose phase-coherence characteristics are incompatible 
with wave-like properties. These are so-called coherent structures, like shocks, vortices and 
convective cells. In fluid turbulence the signature of the presence of coherent structures is 
provided by the existence of non- Gaussian features in the probability density. This is usu- 
ally identified with the velocity-difference probability density function (pdf), traditionally 
adopted for the description of homogeneous turbulence. This explains why in the past the 
treatment of hydrodynamic turbulence was based on stochastic models of various nature. 
These models, which are based on tools borrowed from the study of random dynamical sys- 
tems, typically rely - however - on experimental verification rather than on first principles. 
An example is provided by stochastic models - based on Markovian Fokker-Planck (F-P) 
models of small-scale fluid turbulence recently investigated in the literature by several au- 
thors (including: Naert et al, 1997 2]; Friedrich and Peinke et al, 1999 Luck et al, 
1999 0; Cleve et al, 2000 Q; Ragwitz and Kantz, 2001 6|; Renner et al., 2001, 2002 0,3; 
Hosokawa, 2002 [9]). The validity of phenomenological statistical Markovian Fokker-Planck 
(F-P) models of small-scale fluid turbulence indicate that they are capable of reproducing 
correctly, at least in some approximate sense, key features of the basic phenomenology of 
turbulent flows. Their approach is based on the assumption that the probability density 
associated to the velocity increments should obey a stationary generalized F-P equation. 
Experimental evidence [2j shows reasonable agreement both with the Markovian assump- 
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tion and the F-P approximation, at least in a limited subset of parameter space. However, 
several aspects of the theory need further investigations. In particular, still missing is a 
consistent statistical description following uniquely from the fluid equations. The theory 
should be able, specifically, to recover correctly the structure functions, characteristic for 
the appropriate turbulence regime, but also - possibly - to apply to non-Gaussian probability 
densities. The latter is, in fact, a typical feature suggested by experimental observations, 



performed at sufficiently short scale-lengths in the inertial range 



lOt ] . Based on a recently 



proposed inverse kinetic theory for classical and quantum fluids (Ellero and Tessarotto, 2004- 
2008 Q, Q), a statistical model of small-scale hydrodynamic turbulence is proposed which 
holds for a generic form of the stochastic probability density g(a) . The approach is intended 
to determine the local pdf's (i.e., the local position- velocity joint probability density func- 
tions), both for the stochastic-averaged fluid fields (Z(r,t)) and - unlike Ref. 1] - also their 
stochastic fluctuations 5Z(r,t, a), respectively denoted as (/) and 5f. In particular, it is 
proven that (/) and Sf uniquely determine, by means of suitable velocity-moments, the fluid 
fields (Z(r,t)) and 5Z(r,t, a). Key feature of the approach concerns the construction of the 
statistical evolution equations for (/) and 5f. In particular, it is shown that (/) is generally 
non-Gaussian and obeys an H-theorem. Finally, (/) it is shown to obey - under suitable 
asymptotic assumptions - to an approximate Fokker-Planck kinetic equation which hold in 
principle even in the case on non-stationary, non-isotropic and non-homogenous turbulence. 

II. STOCHASTIC INSE AND STOCHASTIC IKT 

It is convenient first to recall the basic equations for the average and stochastic fluid fields 
and the corresponding initial-boundary value problem. Starting from the incompressible 
Navier-Stokes (NS) equations (INSE) and invoking the stochastic decomposition given above 
the relevant stochastic fluid equations read 

(NV) = (1) 
V • (V) = , (2) 

NV - (NV) = 0, (3) 
V ■ 6V = 0, (4) 
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to be denoted as stochastic incompressible NS equations (stochastic INSE) and N is the 
Navier-Stokes operator iVV =J^V + V • VV + — [Vp — f ] — zA7 2 V. In particular the first 
ones ([I])-© are hereon denoted as stochastic- averaged INSE. These equations are assumed to 
be satisfied pointwise in a set Q x J, being Q an open subset of M 3 and I a finite time interval. 
Assuming that the fluid fields and f are sufficiently smooth, the conditions of isochoricity 
and incompressibility imply the validity of the Poisson equations respectively for (p) and 5p. 
Let us now adopt the statistical approach developed in Ref . [3] , which allows us to cast the 
stochastic INSE problem in terms of a so-called inverse kinetic theory (IKT) [l3j|. This is 
based on the the identification of the fluid fields Z = {V,p} with the moments of a suitable 
probability density /(x, t; Z) defined in the extended phase space T x I [being x the state 
vector x = (r, v) G T spanning the phase-space T = Q x V, with Q the fluid domain and 
V = M 3 the corresponding velocity space]. In particular, it follows that /(x, t; Z) must obey 
an inverse kinetic equation (IKE), which can be identified with the Vlasov-type equation 



L(Z)f = 0. 



(5) 



Here the notation is standard. Thus L is the streaming operator L(Z) = + ^ • {X(Z)} 
and X(Z) = {v, F(Z)} , with F(Z) = F(x, t; f, Z) a suitable vector field (mean field force). 
Subject only to the requirement of suitable smoothness for the fluid fields, several important 
consequences follow |12| . in particular: 



The mean-field force reads: 

F(x, t) f,Z) = — [V-fl - Vpi - f] + \u • VV+^VV • u + z/V 2 V+ 



1 

+ -u 

2 



Po 

^ + V-Vlnp 1 + l(V-Q- [V-nJ-Q) 

at pi L — 1 ' 



1 i ir 

+ — v-n 



Po 



-til 



(6) 



where Q and II are the velocity moments J d 3 vGf for G = u 1 ^ , uu and u=v — V(r,t) i 



is 



the relative velocity. 



{V,p}can be identified in the whole fluid domain Q with the velocity moments G = 
of /(x, t; Z), i.e., respectively, V(r,t) = J <i 3 t>v/(x,t; Z) andp(r,t) =pi(r,t)— P , 



where Po(t) (pseudo-pressure) is an arbitrary strictly positive real function of time 
defined so that the physical realizability condition p(r,t) > is satisfied everywhere in 
Q x I and pi(r,t) = p Q J <iv|u 2 /(x, t; Z) is the kinetic scalar pressure. 
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If f(x,t;Z) is a strictly positive and summable phase-function in T, the Shannon 
entropy functional S(f(t)) = — f r dbcf (x, t; Z) In f (x , t; Z) exists Wt E I and can be 
required to fulfill a constant H-Theorem (see Ref. 14J), i.e., -j^S(f(t)) = 0. 



• Introducing the notation x 2 = v\ h = 2pi/p Q , the pdf 

3/2 

/m(x,*; Z) = V^ ex P i~ x2 } ( 7 ) 

[local Maxwellian kinetic equilibrium) is a particular solution of the inverse kinetic 
equation ([5]) if and only if {V,p} satisfy INSE. 

It is now immediate to obtain an inverse kinetic theory for the previous stochastic 
fluid equations. In fact, let us assume that the operator () is taken at constant r,v 
and t and introduce the stochastic decompositions /(x, t; Z)= (/(x, t; Z)) + 8f(x,t;Z), 
L(Z)=(L{Z))+5L(Z) andF(x,t;/,Z) = (F)+SF, where (L(Z)) = |+v^ + ^. {(F(Z))} 
and respectively 5L(Z) = -j^ ■ {(F(Z))} . Then the following theorem holds: 



Theorem 1 - Stochastic IKT for the stochastic INSE problem 

// /(x, t; Z) is a particular solution of the IKE [Eq. then it follows that: Bl ) (/(x, t; Z)) 
and Sf(x.,t;Z) obey the coupled system of stochastic kinetic equations 

(L(Z)) (f) = - (5L(Z)5f) = (C) , (8) 

(L(Z)) 5f = -5L(Z) {(f) + 5f} + (5L(Z)8f) ; (9) 

B2) A particular solution is provided by (f) = (/m (x,t; Z)) , 5f = dfui^-M Z); B3) Eq. 
can also be written in the form 

(L(Z)) (A/(x, t- Z)) + (AX) /(x, t; (Z)) = - <5AX5A/(x, t; (Z))) , (10) 

where AL = L(Z) -L((Z)), A/(x,t;Z) = /(x,t;Z) - /(x,t; (Z» and L((Z)), /(it; (Z» 
are respectively the stochastic- averaged streaming operator and pdf given in Ref. fil) [see 
Eq.(5)] and there results (C) = — (5AX5A/(x, t; (Z))) ; B4) If Po{t) is defined so that the 
constant H-theorem -j^S(f(t)) = is fulfilled identically in I, it follows that (/(x,£; Z)) satis- 
fies the weak H-theorem §iS({f(t))) > 0, hence both /(x,i; Z) and (/(x,£; Z)) are probability 
densities. 
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PROOF - The proof is immediate. In fact Bi) follows invoking the inverse kinetic Eq.Q 
and the stochastic decompositions for /(x, t; Z), L(F) and F(x, t;f, Z). B2) is mani- 
festly fulfilled since by construction /m(x, t;Z) is a particular solution of the inverse ki- 
netic equation (jSJ). Proposition B 3 ) follows by noting that /(x, t;Z) can be represented 
as f(x,t;Z) = f(x,t;(Z)) + A/(x, t;Z), where f(x.,t;(Z)) is the stochastic-averaged 
pdf solution of the IKE of the form L((Z))/(x, t; (Z)) = , where L({Z)) is a suitable 
streaming operator (see Eq.(5) in Ref.|l|). Finally for B4) we notice that by construction 
f r dx(f{x,t;Z)) = 1. It follows f r dxf(yL,t;Z)\n(f(x,t;Z)) < / r dx/(x,t; Z) ln/(x,i; Z) 
and hence S( fit)) < S((f(t))). Hence, requiring since Po(t) can be determined so that 
-^.S(f(t))=0 14j, this proves the H-theorem. 

An obvious implication (of this result) is the manifest non-Gaussian behavior of the 
stochastic-averaged pdf (f(x,t;Z)). In fact, even if /(x, t; Z) is identified with the local 
Maxwellian distribution /^(x,t; Z), for arbitrary prescribed choices of the stochastic prob- 
ability density g(a) its stochastic average is generally non-Gaussian since (/M(x,t;Z)) 7^ 



III. FOKKER PLANCK APPROXIMATION 



It is interesting to stress that Eqs. ([H]) and are formally similar to the Vlasov equation 
arising in the kinetic theory of quasi-linear and strong turbulence for Vlasov- Poisson plasmas 



and related renormalized kinetic theory [19(, which are known to lead generally 
to a non-Markovian kinetic equation for (/) alone. Nevertheless, the stochastic-averaged 
kinetic equation [i.e., Eq. ([8])] is known to be amenable, under suitable assumptions, to an 
approximate Fokker-Planck kinetic equation advancing in time (/) alone. This is achieved 
by formally constructing a perturbative solution of the equation for the stochastic per- 
turbation Sf. To obtain a convergent perturbative theory, however, this usually requires the 
adoption of a suitable renormalization scheme in order to obtain a consistent kinetic equation 
for (/) . An analogous suggestion is posed by the phenomenological Fokker-Planck models 
of small-scale fluid turbulence recently investigated in the literature. This suggests to seek 
for a possible approximate representation of this type holding for the stochastic-averaged 
kinetic equation (jSj) which should hold for generally non-Gaussian pdf 's (in fact it is obvious 
that generally the average distribution function (/(x, even in the case in which it 
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coincides with (/m(x, t; Z)) results generally non-Gaussian. Let us now first assume that 
(/) = (/a/(x, £; Z)) . In such a case the pseudo-pressure Pq can be defined - consistent with 
Eq. (fT3~l) - so that locally in phase-space the following asymptotic orderings 

J^«i~ (C), (ii) 

X = o(C), (12) 

are satisfied, being ( a dimensionless infinitesimal parameter. In such a case, without 
additional assumptions on the amplitude of the stochastic fields and in the sub-domain of 
velocity space in which (|12|) holds, the "collision operator" (C) in Eq.flT]) can be approximated 
by the generalized Fokker-Planck (F-P) operator of the form 

(13) 

with summation carried out on j, m from to oo or to finite N > 3 and with 
n — j + m > 1, being Cpp a F-P operator and C Jim suitable F-P coefficients. In the 
previous equation /(x, t; (Z)) = /a/(x, t; (Z)) which denotes the local averaged-Maxwellian 
/a/(x, t; (Z)) = - i exp {-X 2 } , with X 2 = g!, v\ h = 2 ( Pl ) / Po and (u) = v- (V(r, t)). 
An analogous result holds also in the case / ^ /Af( x , In particular, the following 

theorem holds: 



(C) ^ C FP 



y- 



gn 



Theorem 2 - Approximate Fokker-Planck IKE 

Let us assume that /(x, £;Z) is a particular solution of the IKE [Eq.$5j)], then provided 
f(x.,t;Z) is a function of (u,p, r,t) which depends slowly both on the relative velocity 
u = v — V and fluid pressure p, in the sense that, in validity of the asymptotic orderings 
( [HP and ( [HP , the local pdf is taken of the form f = f((u,(p,r,t) (smoothness assumption). 
Then in follows that: 

Bl) Eg. (|?3p holds also for a generic pdf (f) ; 

B2) there exists a minimal representation for the F-P operator ( T73j) obtained by re- 
taining only the following F-P coefficients: C^i = — (5F5p) and C it0 = (5F5Ai) with 
i = 2,3,4 where respectively 5A 2 = 5V,5A 3 = -5V6V and 5A 4 = |5V (5V5V) - 
[!<5V5V5V— , such that the stochastic-averaged kinetic equation Eqs. |2P re- 

covers the exact stochastic-averaged fluid equations (CP-([ip. 

PROOF - The proof (of proposition B x ) follows by explicitly evaluating SAf in terms 
of /(x, t; (Z)). This can be achieved formally by introducing a suitable Taylor expansion 

9 



for SAf in the neighborhood of /(x, t; (Z)). This is permitted if / satisfies the smoothness 
assumption (jSj) indicated above, which implies in particular that the orderings ffTTl) . f|T2|) 
must be satisfied. This requirement is manifestly satisfied, in particular, by / = /a/- For the 
general case in which / ^ J'm this permits to estimate SAf in terms of /(x, £; (Z)) by a 
perturbative (Taylor) expansion. In particular, to leading-order in 0(Q, SAf reads SAf = 
~~ "Bv ^PWq) /( x '^i (^)) [1 + °(C)] • Since the n-th order terms in Eq. (jT3l) result by 
construction of order o(( n ), this implies the asymptotic convergence of the expansion for Sf 
in the velocity sub-domain in which the ordering (fl2l) holds. Furthermore, one can prove 
that in order that the stochastic kinetic equations (jHJ) and ([9]) yield a consistent inverse 
kinetic theory for the stochastic-averaged N-S equations (including the corresponding energy 
equation), it is sufficient to retain only a finite number of F-P coefficients (as indicated in 
proposition B2). The proof by explicit calculation of the relevant moment equations of IKE 
in which the truncated F-P collision operator is invoked. 

A basic consequence of THM.2 is that in the velocity-space sub-domain defined by the 
inequality \X\ ~ o(C)> the stochastic-averaged kinetic equation (CQ) can be approximated by 
a time and spatially-dependent F-P equation containing generally both velocity and pressure 
perturbations. In particular, due to the asymptotic ordering ffTTT) . f|T2|) . to leading order 
in o(() pressure perturbations appear only through the explicit contributions carried by Sf 
and Ci,!. Finally, the form of the F-P operator is independent of the specific choice of the 
pdf, provided the above smoothness assumptions are satisfied. 



IV. CONCLUSIONS 



ence has been formulated, based 
Q, Q, which holds for 
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In this paper a statistical model of hydrodynamic turbu 
on the IKT approach earlier developed by Ellero et al. 11, 
a generic form of the stochastic probability density g{a). Basic feature of the new theory 
is that it satisfies exactly the full set of stochastic fluid equations while permitting, at the 
same time, the construction of the stochastic pdf which - in difference with Ref.|l| - advances 
in time the full set of stochastic fluid fields. Unlike customary statistical approaches, this 
is identified with the local position-velocity joint probability density function, rather than 
the two-point correlation function (velocity-difference pdf). The present theory displays 
several interesting features. In particular, the stochastic-averaged pdf has been shown to 
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17l | , for which a renormalized kinetic theory can be in principle applied 



satisfy an H-theorem, assuring its strict positivity. The corresponding kinetic equation is 
formally similar to the Vlasov equation arising in the strong turbulence theory of Vlasov- 
3 oisson plasmas 

1^ |. However, unlike plasmas, the stochastic IKT here considered admits exact local 
kinetic equilibria for (/(x, t; Z)), which are expressed by the stochastic-averaged Maxwellian 
distribution t; Z)). More generally, for pdf's which are suitably smooth, in the sense 

of Thm.2, the stochastic-averaged kinetic equation can be approximated in terms of Fokker- 
Planck kinetic equation. This result suggests a possible new interesting viewpoint for the 
investigation of turbulence theory in neutral fluids, which includes, in particular, the analysis 
of previous F-P models of turbulence [2|. This topic will be the subject of a forthcoming 
investigation. 
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